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Abstract 

We compute the leading non-analytic contributions of the form m q log m q to 
matrix elements of twist-2 operators in the nucleon and pion using effective 
. field theory. Previously omitted one-loop contributions that are related to 

tree- level matrix elements by chiral symmetry are included. 
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I. INTRODUCTION 



Deep-Inelastic-Scattering (DIS) from nucleon targets has provided a wealth of informa- 
tion about the nature of strong interactions and the structure of the nucleon. From the 
initial discovery of partons in the 1970's, DIS continues to provide ever more precise mea- 
surements of the parton-distribution functions (PDFs). It is highly desirable to make a direct 
connection between the experimental data and the now well-established theory of strong in- 
teractions, QCD. A rigorous and model-independent comparison will be accomplished by 
performing high-statistics unquenched or partially-quenched |J lattice-QCD calculations [§. 
Presently, lattice calculations cannot be performed with the physical values of the light quark 
masses, m q , (m u ~ 5 MeV, ~ 10 MeV) and extrapolations from the lattice masses, that 
produce a pion of mass mj^' ~ 500 MeV, to the physical values must be performed. Of 
course, such extrapolations require knowledge of the m^-dependence of the matrix element of 
interest. Recently hadronic models, such as the Cloudy Bag Model (CBM), have been used 
to motivate explicit forms for the m g -dependence of forward matrix elements of the non- 
singlet twist-2 operators that contribute to DIS from the nucleon @,f§]. In addition, these 
models have been used to connect lattice calculations to other properties of the nucleon, 
such as electromagnetic form factors || . 

It is well established that one can determine the m 5 -dependence of hadronic observables 
by performing a systematic expansion about the chiral limit 0-^. In fact, extensive work has 
been accomplished in understanding the properties and interactions of the low-lying mesons 
and baryons in both two-flavour and three-flavour QCD, such as the magnetic moments (JTO 



the electric form factors [[11]], the axial matrix elements and the polarizabilities of the 
nucleons and other octet baryons to name just a few, and the analogous quantities in 
hadrons containing heavy quarks |L4] . In addition, there has been an extensive effort during 



the last decade to include multi- nucleon systems in this framework ]T5|. In this work, we 



include twist-2 operators in the chiral lagrangian, and compute the leading non-analytic 
contributions of the form m q log m q to their matrix elements. 

The pion fields are introduced into the low-energy effective field theory (EFT), chiral 
perturbation theory (%PT), through the E-field, 

with / = 132 MeV. E = £ 2 transforms as E -> LEi?t = l^W U£R? under SU{2) L ®SU{2) R 
chiral transformations. In order to construct an EFT with well-defined power-counting, 
the nucleons are treated as heavy fields in the Heavy-Baryon formulation of Jenkins and 
Manohar , and transform as N v — > U N v (the subscript v denotes the four- velocity of the 
nucleon) under chiral transformations (for reviews see Ref . |m ) . Below the chiral symmetry 
breaking scale A x , S-matrix elements can be expanded in derivatives and in m q . The naive 
size of the matrix element of an operator with n\ creation operators for heavy nucleons, n 2 
annihilation operators for heavy nucleons, n% derivatives, n 4 light quark mass matrices, n 5 
powers of the nucleon four-velocity v, n 6 E-field operators and n 7 E^-field operators, is 




Ax, 



Er(Et) 7 , (2) 
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which should be considered merely as a guide. The strong interactions between pions and 
nucleons at leading order in the chiral expansion arise from a lagrange density of the form 
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f 

£ = ^-Tr 



} + A Tr [ m E t + h.c 



"q <- 



+ N v iv ■ DN V + 2g A N v S ■ A N v , (3) 

where = + V fJ- is the chiral- covariant derivative with V M = \ {^d^ + C^/iC) the pion 
vector field, g& = 1.25 is the axial- vector coupling constant, is the covariant spin-operator 
defined in Ref. A^ = | f^9^ — C^/iCJ is the axial- vector pion field, and A is a parameter 
that provides the leading order contribution to the pion mass. 

II. NON-SINGLET, TWIST-2 OPERATORS 

In this section we will focus on forward matrix elements of the non-singlet twist-2 oper- 
ators 

... Mn = ^1 r a 70*i (i B w ) ... (i q - traces , (4) 

where the {...} denotes symmetrization on all Lorentz indices. Once the transformation 
properties of 0<$£ ^ under SC/(2) L ® SU(2) R have been established, the operators in 
the pion-nucleon EFT that reproduce matrix elements of O^'^ Mn can be constructed, 
and further, the EFT can be used to perform a model-independent calculation of the low- 
momentum contributions to its matrix element. In order to determine its transformation 
properties it is convenient to write it in terms of left-handed and right-handed quark fields, 

Wn),a _ /*}(n),a . £j(n),a /p-x 



with 



°l^2 ...n = — { Ql T l 7{ w (iD^ 2 ) ... (i D^ } ) q L - traces 



1 - _., 



°jTmia« .../*„ = ^[ ^ T R 7{mi A«a ) - (* ) - traces , (6) 



from which it is clear that 0^ a 2 ^ n transforms as (3, 1) © (1, 3) under SU(2) L ® SU(2) R . 
For bookkeeping purposes we have introduced the flavour matrices r£ and r R that are taken 
to transform as r£ -> Lr£Lt and r£ -> i?r^i?t under S77(2) L ® SU(2) R . 

If we are interested in DIS from pions, or in pion loop contributions to DIS from the 
nucleon we require the matrix element of ^ n in the pion. The power-counting rules of 

eq. ([|) dictate that such matrix elements will be dominated by operators involving the least 
number of derivatives and insertions of m q , but can have an arbitrary number of insertions 
of the S and fields. Operators that can contribute are of the form 

..... "> a(n) «" X Tr [SW W ...^E + Sr^ w ...^St] - traces 

= a( n h (i) n ie aa ? vr a ^ 1 ...^ n 7r /3 - traces + (vr 4 ) . (7) 
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With the exception of n = 1, the coefficients are unknown and must be determined 
elsewhere. is the isovector charge operator from which we deduce that aS l > = +1. In 
addition to the operators of eq. ([?]), there are also operators of the form 




O l >[):: .. <i u " (/)"—( — ) Tr [E ^ fc E^ fc ...^ p Et ^..."^ n} E 

traces , (8) 

that must be considered. However, these operators do not contribute to single-pion forward 
matrix elements at tree-level, and only give interactions between three or more pions. In 
addition, the symmetry of the Lorentz indices and the tracelessness of O&ll" means that 
these operators do not contribute to single-pion forward matrix elements even at one-loop 
level. Therefore the diagrams shown in fig. [I] give the leading non-analytic contributions 



; at 



(a) (b) 




(c) 

FIG. 1. The pion loop diagrams that give the leading non-analytic contributions to the matrix 
element of C^")]" ...^„ between single-pion states. The crossed circle denotes an insertion of an 
operator from eq. (Qj, arising directly from the twist-2 operator. The smaller solid circle denotes 
an insertion of a leading order strong-interaction vertex from eq. Diagrams (a) and (b) are 
vertex corrections while diagram (c) denotes wavefunction renormalization. 



to the matrix element of 0^'" Mn between single-pion states. After a straightforward 
calculation, one finds that the forward matrix element of 0,\ n };° ,, between an initial pion 
with isospin index a and momentum q^, and a final pion with isospin index (3 vanishes for 
n-even, and is 

e afia q^-q^ - traces , (9) 

for n-odd, where the ellipses denote terms that are analytic in m q , or are higher order in 
the chiral expansion. The factor of 1 — 5 nl in the sub-leading contribution ensures that the 
n — 1 isospin-charge is not renormalized at loop-level. 
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The leading order operator contributing to the matrix element of ^ n between 

nucleon states is 

...^ - AW v^.-.v^ N v r* + N v - traces , (10) 

where operators involving more derivatives or insertions of m q are suppressed by powers of 
A x , and we have defined rj?± to be 

r! ± = ~ ( ± ^ ) . (11) 

The coefficients must be determined elsewhere, except for = +1 which corresponds 
to the isospin charge operator. In addition to the tree-level contribution to the nucleon 
forward matrix element, there are vertices involving an even number of pion fields. The 
diagrams shown in fig. |^ give the leading non-analytic corrections to the forward matrix 
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FIG. 2. The pion loop diagrams that give the leading non-analytic contributions to the matrix 
element of ...^ n between single-nucleon states. The crossed circle denotes an insertion of an 

operator from eq. (f^) or eq. (]T^) ; arising directly from the twist- 2 operator. The smaller solid circle 
denotes an insertion of the strong two-pion-nucleon interaction from the nucleon kinetic energy term 
in eq. (fjj, while the square denotes an insertion of the axial-vector interaction oc gA- Diagrams 
(a)-(d) are vertex corrections while diagram (e) denotes nucleon wavefunction renormalization. 



element of O^'^ ^ between single-nucleon states. A straightforward computation gives 



M=AW v^.-.v^ U v r a U v 



1 - 3^ + 1 



1-5 



rat 



l% 2 f 2 



-m\ log — y 



A 2 

X 



(12) 
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where U v is the nucleon spinor, and we have only shown the non-analytic part of the sub- 
leading contribution. The non-analytic contributions to the n = 1 matrix element vanish as 



the nucleon isospin charge is not renormalized. Our result in eq. (|T2| ) differs from previous 



computations @,|] primarily due to our inclusion of the two-pion-nucleon interaction asso- 
ciated with the tree-level matrix element of O.Y?};" „ , but also in the numerical coefficient 
of the g\ contribution. In Ref. || only the contribution from diagram (b) of fig. |2| is com- 

71 + 1 

puted and correctly found to scale as m q 2 \ogm q . However, this diagram is only part of the 
complete result for n = 1 (required by charge conservation) and is sub-dominant for n > 1. 
In Ref. [f|], this was corrected somewhat by the appearance of m q \ogm q contributions for 
all values of n, however, the qa independent contributions were omitted. 

It is worth keeping in mind the relative size of non-analytic terms from loop-diagrams 
compared with the analytic contributions from both loop-diagrams and local counterterms. 
The complete set of local operators with a single insertion of m q that contribute to the 
matrix element of 0[! l }] a ' ,, is 



£41/^2 



0%>£ ..., n - ( 6i(A*) N {t« + , X+ } N + bM N [r« + , x-\ N + b 3 (») Tr [ X+ ] N r« + N 
+ h(fi) N {r*_,x-} N + bM N [t?_,x+] N ) V/J2 ...^ n , (13) 

where £± = ~ ^m,(t ±£mj£^. In general, the constants 6«(/i) are renormalization scale 
dependent, and must be determined experimentally. Only the operators with coefficients b\ 
and 63 contribute to forward matrix elements in the nucleon. The operators with coefficients 
62 and 65 have at least one additional pion associated with them, while 64 has at least two 
additional pions associated with it. The m g -dependence of the matrix elements of O™" Mn 
shown in eqs. @, ([12]) and (|D|) has the form 

M ~ a + (3m q log (rn q /fi) + ^(fi)m q + ... , (14) 

where the /i-dependence of 7^) is precisely equal and opposite that of the /5-term, leaving 
an expression that is explicitly \x- independent. While the contributions from the /3-term 
formally dominate the sub- leading contributions in the chiral limit when fi = A x , there are 
well-known examples where such terms are numerically smaller than the sub-leading analytic 
contributions, terms analogous to the 7-term, for the physical values of m q . The pion-charge 
radius (r£) is such an example, where for the physical values of the pion mass (and kaon 
mass in 577(3)), the contribution from the ag(fi) counterterm (evaluated at // = A x ) is twice 
that of the non-analytic loop contributions of the form logm q . Therefore, while the terms 
we have computed in eq. ( |T2] ) are the formally dominant sub-leading contributions in the 
chiral limit they may not dominate the sub-leading contribution for physical values of m q 
due to the terms shown in eq. flT3j). 

At relatively low momentum scales, ~ 300 MeV, there can be large contributions from 
loop diagrams involving the A's. The formal construction and phenomenology of an EFT 
with dynamical A's (or any resonance) has been studied extensively [|7|- [T^J16[| . If all diagrams 
with the A-resonance as an intermediate state are included then one can consistently take 
/i ~ A x , and capture the dominant infrared behavior of the theory J7LH|. Matrix elements of 

are described, at leading order, by 
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Wn),a 

— /J-n 



C (n) v^v^.-.v^ A" t£ + A ai „ - traces 
1 



+ D {n) —v^v^.-.v^ A^_ liV rg + A Mn}> „ - traces 



(15) 



where C^ 1 ' = —3 by normalization of the isospin charge operator, and = simply 
because of the number of available Lorentz indices. These operators will contribute to 
matrix elements between nucleon states through loop diagrams. The leading order strong 
interactions between the iV's, A's and 7r's are described by a lagrange density of the form 



9na 



A a ' ijk A l ak N, e u + h.c. 



(16) 



where the coupling is g^A ~ 1-8, from the observed width for A — > Nn ||. The loop 
diagrams shown in fig. ^| gives rise to a forward matrix element between single-nucleon 
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FIG. 3. One-loop diagrams with A intermediate states that contribute to the matrix elements of 
0^™)i2 —tin between single-nucleon states. The thick solid line inside the loop denotes a A propaga- 
tor, while the dashed line denotes a pion propagator. The crossed circle denotes an insertion of an 
operator from eq. ^ or eq. (It ), arising directly from the twist-2 operator, and the square denotes 
an insertion of the strong NAtt interaction oc gNA- Diagrams (a) and (b) are vertex corrections 
while diagram (c) denotes nucleon wavefunction renormalization. 



states of 

M -- 

where AM 
Ji(A,m) : 
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4tt 2 / 2 



J x (AM, uin) U v r a U v 
- traces ) 



A (n) , 5 ( n) _ _5_£(„) , 2 1 

9 27 3 



Ma — Mtv is the A-N mass difference. The function Ji is 

,2\ / / 



m 2 - 2A 2 ) log ( ) + 2AVA 2 - m? log 



v/A 2 



m 2 + ie 



A 2 

. x 



A + VA 2 - m 2 + ie 



(17) 



(18) 



where we have renormalized at the scale \l — A x . In the limit that A — > 0, J\ contains a 
chiral logarithm, as is clear from eq. (fOf), while in the limit of large A only terms analytic 
in m survive, as required by the decoupling of the A ||. The reason that such contributions 
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must be included in order to sensibly renormalize the theory at /i = A x is that the scale for 
the m^-dependence from these diagrams is set by AM and not by A x , and therefore a naive 
estimate of the size of counterterms in the theory without the explicit A-fields is set by AM 
and not A x . Without resorting to hadronic models, one is unable to make statements about 
(ji n ) or £)( n ) ) anc i they must be determined elsewhere. It is interesting to note that NA 
transition operators induced by 0^' 2 Mn involve either a spin-operator S 1 * or additional 
pion fields. Neither type of operator contributes to nucleon matrix elements at one-loop 
level. Our result in eq. (|17D disagrees with the result of Ref. 0, as they computed only the 
contribution from fig. ^| (b), where Oj$s£ is inserted into the pion propagator. 



III. SINGLET, TWIST-2 OPERATORS 



In this section we consider matrix elements of singlet twist-2 operators, ^ s 'Ojj^ 2 ^ with 
j = q,g for the quark and gluonic operators, of the form 



{S) Ofl ...Mn = ^ V 7{ w (i D» 2 ) ... (i 5 Mri} ) q - traces 

{S) °& ..... = h G k»i (* A*) - (i D^) G% - traces . (19) 



and it is clear that the {s) Oj}$ 2 ^ n transform as (1, 1) under SU(2) L g> SU(2) R . 

The matrix element of ^Ojj^ ^ in the pion will be described at leading order by a 
lagrange density of the form 



c2 ' • \ n ~ 1 



(S) Ol^l - af ] (if f -{^\ Tr Et^...^S + E^...^ 



traces 



/ 1 \ n_1 

= afh (iY f F J ^d^.-.d^ - traces + (V) , (20) 

where the a,j are coefficients that must be determined elsewhere, and depend upon the par- 
ticular singlet operator under consideration. A calculation of one-loop diagrams analogous 
to those shown in fig. |T] give rise to a matrix element between pions with isospin indices a 
and P, at next-to-leading order, of 



-y \ n—l 



M j = 4a? ) l — \ 5^ q^.-.q^ - traces , (21) 

for n-even, while the matrix element for n-odd vanishes. There are no non-analytic correc- 
tions to this matrix element for any n at next-to-leading order. 

In contrast to eq. (|T0"|), the leading-order operator contributing to the matrix element of 
a singlet operator ^'OV^L ...„ n between nucleon states is 

{S) OS 2 ..., n - v^.-.v^ N v N v - traces , (22) 

where operators involving more derivatives or insertions of m q are suppressed by powers 
of A x . Since is the baryon number operator and ^O^j} 2 vanishes for n < 2, we 
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have that A q V = +3 and = 0. Only some of the one-loop diagrams in fig. |2] contribute 
to matrix elements of the singlet operators. Diagrams (c) and (d) of fig. |2| are absent 
while diagram (b) can only contribute at higher orders. Further, the contribution from 
the vertex diagram, diagram (a), is exactly canceled by the contribution from wavefunction 
renormalization, diagram (e). Therefore, the singlet matrix elements in the nucleon do 
not receive any non-analytic corrections of the form m q log m q from nucleon intermediate 
states [[ However, there are contributions from A intermediate states. The leading order 
matrix elements involving the A are described by the operators 



(•SVoiO) _> n( n ) v v v A a A — traces 

+ Df ] —(v { ^v^...v^ n _ 2 A^ uV A Mn})U - traces , (23) 



with = —3, C^p = 0, = 0, and Djp = 0. These operators contribute through 
loop-diagrams shown in fig. |3|. Diagram (b) of fig. [3] contributes only at higher orders in the 
chiral expansion, and we find a contribution to the nucleon matrix element of 



9na 
4tt 2 / 2 

( w M i^ 2 ...v - traces ) . (24) 



As required, corrections to the matrix element of the baryon number operator, ^'O^ 1 ', 
vanish. 



IV. CONCLUSIONS 

We have computed the leading non-analytic contributions of the form m q log m q to matrix 
elements of twist-2 operators that arise in deep-inelastic scattering. A previously omitted 
contribution to the matrix elements of non-singlet operators that is independent of qa and 
related by chiral symmetry to the tree-level vertex is identified. Our results will aid in the 
extrapolation of unquenched lattice calculations of single-nucleon matrix elements of twist-2 
operators from the quark masses used on the lattice to their physical values. 

The work we have presented here can be straightforwardly extended to off-forward matrix 
elements of the twist-2 operators, i.e. those in which there is a non-zero momentum transfer 
to the hadronic system from the twist-2 operator. Such matrix elements have received 
significant amount of attention during the past few years, as one can define off-forward 
parton distributions as a simple extension of the parton model (for an overview see Ref. [|TT|1). 
In addition, deeply-virtual-Compton-scattering (DVCS) has been extensively explored as 
a means to measure such distributions. The effective field theory construction we have 
employed in this work will allow for a description of these matrix elements in the low- 
momentum regime. 



The loop corrections are the same as those contributing to the nucleon mass, for which there are 
no terms of the form m q logm q . 
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